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Evolution of a spinor condensate: coherent dynamics, dephasing and revivals 
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We present measurements and a theoretical model for the interplay of spin dependent interac- 
tions and external magnetic fields in atomic spinor condensates. We highlight general features like 
quadratic Zeeman dephasing and its influence on coherent spin mixing processes by focusing on a 
specific coherent superposition state in a F — 1 87 Rb Bose-Einstein condensate. In particular, we 
observe the transition from coherent spinor oscillations to thermal equilibration. 
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Multi-component Bose-Einstein condensates with spin 
degree of freedom, the so called spinor condensates, are 
experiencing rapidly growing attention. These ultra-cold 
quantum gas systems are interesting in several respects, 
e.g. they show intriguing static and dynamic magnetic 
properties Q, 0, El , they represent a well controlled ther- 
modynamic model system 0, 0] and they promise unique 
insights into mesoscopic multi-component entanglement 
and decoherence processes @. 

Pioneering experiments on F = 1 23 Na spinor con- 
densates, found to be anti-ferromagnetic, and quasi- 
spin-1/2 systems in 87 Rb have shown fascinating demix- 
ing dynamics, metastability and domain formation pro- 
cesses 0, 0- Recently the spinor systems under study 
have been extended to 87 Rb, which behaves ferromag- 
netic in the F = 1 state 0, and anti-ferromagnetic 
in F = 2 Q. In particular, these experiments also 
showed spinor oscillations and shifted the research in- 
terest towards dynamic spin conversion effects. Recent 
studies have demonstrated a magnetic field dep endence 
of the oscillation amplitude and frequency 0, 0> E3 j 
which was also theoretically assigned to an interesting 
resonance phenomenon [l3j . The coherent evolution in 
these experiments at short timescales is opposed to the 
incoherent thermodynamic behavior at long timescales, 
e.g. decoherence driven cooling Q, constant temperature 
BEC |f| or temperature driven magnetization 

The coherent and thermodynamic regime are concep- 
tually very different, as the first relies on the assump- 
tion of each atom in the ensemble being in the same 
superposition state, while there are presumably differ- 
ent and independent ensembles for each spin state in 
the second regime, i.e. here an atom is cither in states 
\F = l,mp = ±1) or \F = l,m_p = 0). It is still under 
discussion, how and at which point the crossover from the 
"pure state" ensemble to a "mixed state" arises and how 
this is connected to decoherence and dephasing mecha- 
nisms in spinor systems. This question can be generalized 
to the understanding of decoherence in arbitrary multi- 
component systems and its dependence on the number of 
constituents, which might have important consequences 
for quantum information theory. So far studies on the 



relative phases of different components have only been 
performed in quasi-spin 1/2 systems 0, 0, based 
on the preparation of mixtures with well defined phase, 
while investigations on spinor systems were essentially 
restricted to population-based state preparation without 
phase control. 

In this paper we discuss dephasing and decoherence 
effects in spinor quantum systems. In particular wc 
present measurements on the influence of finite tem- 
perature and external magnetic fields on the evolution 
of a F = 1 87 Rb spinor condensate prepared in the 
C f = (C+,Co,C-) T = (-l,-*V2,l) T /2 state. In addi- 
tion we develop an analytic description for the mean-field 
evolution and emphasize the importance of the quadratic 
Zeeman effect, which leads to periodic dephasing and 
rephasing of the spin expectation value. We theoreti- 
cally predict and experimentally find dephasing induced 
spin dynamics with magnetic field dependent period and 
amplitude, which clearly shows the coherence of spin mix- 
ing in spinor condensates. An analysis of the data again 
independently confirms the F = 1 state of 87 Rb to be 
ferromagnetic. 



I. THEORY 

In order to simplify the discussion and concentrate on 
the spinor physics, we will restrict ourselves to the single 
mode approximation (SMA), which effectively factorizes 
the F = 1 spinor order parameter as 
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(1) 



where the volume density n(r) and spinor £ are indi- 
vidually normalized to J d 3 rn(r) — N and ||C|| 2 = h 
respectively. This assumption is justified, if the spin dy- 
namics is slow compared to the timescale of motion in 
the trap and if the spin dependent interaction is much 
smaller than the spin independent one. 
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In SMA,the mean-field equations of motion for the 
spinor are [Ella] 

id t C = (H z + H mf )C, (2a) 
H z = H LZ + H QZ = -pF z -q(t-F 2 ), (2b) 

ffmf = .92 (n) ^ ^aCC^a, (2c) 
a— x.y,z 

where p oc B and q oc B 2 are the coefficients of the lin- 
ear and quadratic Zeeman effect respectively. For 87 Rb, 
F=l, the Breit-Rabi formula allows to calculate them as 
p « —^i B B/2h and q w p 2 /whfs with the 87 Rb ground 
state hypcrfinc splitting Whfs ~ 27r ■ 6.835 GHz. The spin- 
dependent mean-field interaction is proportional to g% = 
4irh/m x (a F=2 - a F =o)/3 Q and (n) = J d 3 rn 2 (r)/N. 
We will identify the 3x3 matrices {F x , F y , F z } as the 
standard representation of the F = 1 angular momentum 
algebra 01 ■ I* 1 principle, the mean- field Eqs. could 
be derived also from a Hamilton- Jacobi theory. The con- 
servation of spin and energy follow naturally from ro- 
tational and time-translation symmetries. Due to these 
constraints, it is indeed possible to find integrablc solu- 
tions to the full nonlinear equations of motion. In case 
of the initial state the time-dependent spinor am- 
plitudes can be written in terms of Jacobi elliptic func- 
tions [Isj j 

cn fc (f)dn fc (fy 
fc S 4(f) 
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where s = exp(— i(g2(n) — q)t/2) is a dynamic global 
phase and k = g2 (n) / q is the ratio of interaction energy 
to quadratic Zeeman effect. For the directly measurable 
populations, the solution simplifies to 

ICoWl 2 = (l-ksnl(qt))/2, (4a) 
\C±(t)\ 2 = (l + ksnl(qt))/4. (4b) 

For small |fc| <C 1, the Jacobi elliptic functions can 
be approximated by ordinary trigonometric ones, i.e., 
snfe(x) « sin(x), cn^(x) « cos(cc), drik(x) « 1. Thus, 
the populations of Eqs. Q oscillate with an amplitude 
given by k, and a period which is n/q for small |fc|. 
This solution offers a clear way to determine the mag- 
netic properties of the system as k is negative/positive 
for ferromagnetic/anti-ferromagnetic systems, resulting 
in an increase/decrease of |(o| 2 at the beginning of the 
oscillation (if q > 0). 

At this point, we would like to emphasize the intri- 
cate connection between quadratic Zeeman dephasing 
and spin mixing dynamics: for zero quadratic Zeeman 
shift q = 0, the state Cf docs not evolve according to 
Eq. |4"|. except for a trivial Larmor precession around 
the magnetic field vector, which can be compensated by 



choosing a rotating frame of reference. The intuitive rea- 
son for this lies in the rotational symmetry of the spin 
interactions, leading to a conservation of magnetization. 
In a frame rotated by 90°, the above state is just the 
fully stretched state £ = (0,0, 1), which is an eigenstate 
of H m { and consequently does not evolve. In fact, any 
rotated stretched state belongs to this class of stationary 
states of H m { [23l |. 

The quadratic Zeeman effect however destroys the fully 
stretched state, removing it from the class of stationary 
states, and thus allowing interaction driven spin dynam- 
ics. Two regimes can be distinguished in the resulting 
population oscillation: the interaction dominated regime 
at low magnetic fields, where the oscillation period is es- 
sentially given by the intrinsic spin interactions, and the 
quadratic Zeeman regime, where it is given by the respec- 
tive beat period 7r/q. In both limits k — > and q — > 0, 
the amplitude of population oscillations goes to zero. In- 
terestingly the crossover between the two regimes shows 
a resonance-like structure, as also recently predicted for 
a more general case |l3j . 



II. EXPERIMENT 

In our experimental setup described in we initially 
prepare degenerate 87 Rb ensembles containing up to 10 6 
atoms in the \F = 1, tuf = —1) state in an optical po- 
tential with trapping frequencies uj x : uj y : uj z w 1 : 7 : 40 
and uj x — 2ir ■ 16 . . . 25 s _1 . Using a radio frequency ir/2- 
pulse of typically 40/xs duration, we evolve the spinor 
into the superposition state £(0) = C?-- In rea l space 
this corresponds to a 90° rotation of the classical spin 



(3b) expectation value T = C^C from T = — e z to J 7 = 



A. Phase evolution 

Using the classical picture it is evident, that an- 
other 7r/2 radio frequency pulse will rotate the spin vec- 
tor upward to = +e z or equivalently the spinor to 
C = (1,0, 0) T , transfe rring all the population into the 
\F = l,mp = +1) state. The addition of a delay time 
in between the two pulses results in the spin-1 analo- 
gon of the famous Ramsey experiment for spin- 1/2 par- 
ticles |l9j . which is sensitive to the quantum mechanical 
phase evolution of the system. The last 7r/2-pulse creates 
superpositions of the spinor components and thus trans- 
forms the phase information into population information 
accessible to detection. Thus, the state of the Ram- 
sey spinor is Cr(^) = exp (— iir/2 F x ) exp {iujRpt F z ) Q{t) 
with Q(t) from Eqs. El 

In the following we will analyze the time-dependent 
population distribution in such a Ramsey experiment 
in order to gain insights into dephasing and decoher- 
ence mechanisms. The corresponding analytical results 
in Fig. ^ clearly show a magnetization oscillation with 
the detuning between the Larmor precession frequency 
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FIG. 1: Theoretical prediction for the outcome of a Ram- 
sey experiment with an interacting spin-1 spinor conden- 
sate, plotting the normalized magnetization in the direction 
of the quantization axis |Cr+| 2 — IC-R-| 2 an d the itif = 0- 
projection |Cro| 2 as a function of pulse delay time. The pa- 
rameters were chosen to be \p\ — 2tt ■ 772 kHz (corresponding 
to \q\ = 2vr • 87.2 Hz), k = -0.03 and lu rf = 2tt • 773 kHz for 
the radio frequency field driving the 7r/2-pulses. 



and the radio frequency for driving the 7r/2-pulses, as 
one would also expect from spin-1/2 particles. 

There is however an interesting "beat" feature with 
twice the quadratic Zeeman frequency adding new 
physics to the spin-1 system. For the given parameters 
there are only very small interaction effects barely visi- 
ble as a modulation of the | Co | 2 -oscillation envelope and 
practically undetectable in the magnetization. 

In order to draw a simple picture of the correspond- 
ing physics, it is thus sufficient to neglect interactions 
and to concentrate on the evolution of a single spin- 
1 particle in an external magnetic field aligned along 
the z-axis. In this picture the linear Zeeman shift of 
Eq. (|2b|) adds a time dependent phase of equal magni- 
tude but opposite in sign to the £+ and £~ components: 
Clz(*) = expHt J ff LZ )Cf= (-e- i f t ,_-iV2,e^) T /2. 
In real space this corresponds to a spin vector J- = 
(sin (pt), cos {pt) 1 0) T rotating around the magnetic field 
axis, i.e. the Larmor precession of a magnetic moment, 
just as one would classically expect. 

A different "non-classical" behaviour arises from the 
quadratic Zeeman effect, which corresponds to adding a 
time dependent phase to the Co component: Cqz(^) = 
exp(-itif QZ )<f= (-l,-i\/2e I< ? t ,l) T /2. The resulting 
spin expectation vector T = (0, cos (qt), 0) T now peri- 
odically disappears and reappears, always pointing along 
the original direction. In a Ramsey experiment, the com- 
bined linear and quadratic Zeeman effect show up as a 
periodic dephasing and rephasing of the rotating spin 
vector just as the analytic result shown in Fig. ^ which 
can also be interpreted as a beat note between the two 
slightly detuned transitions mp = <-> mp = ±1. 

Fig.|Ushows the experimental results of such a Ramsey 
sequence consisting of a 7r/2- pulse for state preparation, a 
variable delay time and another 7r/2-pulsc before Stern- 
Gerlach separation and imaging of the individual spin 
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FIG. 2: Magnetization amplitude and trif = population 
after a Ramsey sequence as a function of the delay time be- 
tween Ramsey pulses. The dots (blue) show the experimen- 
tal results for the condensate fraction, while the circles (red) 
represent the normal component. The lines are phenomeno- 
logical envelope fits with exponential decay to guide the eye. 
The magnetic field parameters are the same as in Fig. 



components. In order to emphasize the effect of quadratic 
Zeeman dephasing instead of the trivial Larmor preces- 
sion, we plot the magnitude of the magnetization vector 
instead of the magnetization itself and concentrate on the 
envelope of the corresponding data points (24|. In order 
to achieve the parameters used for Fig. 2] an offset field 
of f=a 1.1 G was applied, corresponding to a Larmor pre- 
cession frequency of \p\ = 27r • 772 kHz and a quadratic 
Zeeman beat frequency \q\ — 2ir ■ 87.2 Hz correspond- 
ing to a beat period of ~ 5.7 ms, i.e. operating in the 
quadratic Zeeman regime. In this case already the single 
particle picture discussed above clearly reproduces our 
experimental data: The magnetization is modulated in 
amplitude by the quadratic Zeeman effect with collapses 
and revivals corresponding to a beat period of ir/q. 

The timescale of the graph was chosen to capture 
the evolution from the coherent to the thermodynamic 
regime, which is reflected by the decreasing envelopes of 
magnetization and mp = population. In order to iden- 
tify further dephasing and decoherence timescales we will 
now concentrate on the evolution of the mp — popu- 
lation, which is not visibly affected by quadratic Zeeman 
dephasing. The oscillation of both normal and condensed 
mp = population is damped and initially tends to- 
wards its mean value of 0.25. This damping could be 
explained by random or spatial dephasing e.g. due to 
magnetic field gradients [20j, or by breakup into small 
dynamical spin domains due to a dynamical instabil- 
ity [2l|. On this timescale, dephasing should at least in 
principle be reversible. In contrast, the drifting apart of 
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normal and condensed mp — population at ps 15 ms 
marks the crossover to the thermodynamic regime on 
longer timescales. In this regime, the system tends to 
restore thermal equilibrium meaning in particular equal 
population (i.e. 0.33) of all spin components in the ther- 
mal cloud 0, 0|- The main process here is redistribu- 
tion of atoms between condensed fraction and thermal 
cloud within each spin component, leading to the diverg- 
ing trend in itif = populations observed. 



B. Population evolution 

In order to observe the effect of spin interaction be- 
yond the single particle level, we investigated the re- 
sulting population evolution of the intermediate Ramsey 
state Qil directly, i.e. in a sequence of a 7r/2-pulse and a 
variable delay before detection. In contrast to the above 
Ramsey experiment this sequence is not sensitive to the 
relative phases of the spinor components but gives com- 
plementary information on the interaction driven popu- 
lation evolution. 

The quadratic Zeeman dominated regime at high mag- 
netic fields is particularly interesting for the investigation 
of coherence in interaction driven spin dynamics. Our 
theoretical results Eq. predict that the quadratic Zee- 
man phase determines the direction of spin mixing dy- 
namics, demonstrating the phase sensitivity of the pro- 
cess. 

The results shown in Fig. [3] clearly confirm the co- 
herence o f sp in dynamics in spinor Bose-Einstein con- 
densates I23. Furthermore the positive initial slope of 
the population evolution of the mj? = state confirms 
the ferromagnetic behavior of 87 Rb in the F=l hypcrfinc 
manifold [113. 

Comparison with theory (Eq. 0} shows an excellent 
agreement during the first 10 ms of evolution, in which a 
coherent behavior is expected considering the results of 
the Ramsey experiment. The crossover to the thermo- 
dynamic regime at later times goes along with damping 
of the population oscillations and a transfer of particles 
between the condensate fraction and normal component. 
The thermodynamic nature of this process also shows up 
as a temperature dependence with longer decoherence 
times for the lower temperature data. It is important to 
note, that the damping of the oscillations observed in this 
case is fundamentally different from the Ramsey experi- 
ment. The Ramsey experiment is sensitive to both spa- 
tial dephasing and thermal decoherence, while the " single 
7r/2-pulse experiment" is only affected by thermal deco- 
herence. The relatively long lasting mp = oscillations 
in this experiment thus corroborate our interpretation 
that the faster decay in the Ramsey experiment is not 
of thermal origin, but probably due to magnetic field de- 
phasing or more complex effects pi) . 

The amplitude of the population oscillations is propor- 
tional to the spin dependent interactions via the param- 
eter k in Eq. (0J. In order to compare the experimcn- 
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FIG. 3: Evolution of the rxiF = population in the con- 
densate fraction for ensembles containing w 10 6 atoms with 
25 ... 30 % condensate fraction (upper figure) and w 4 ■ 10 5 
atoms with 75 ... 80 % condensate fraction (lower figure), pre- 
pared in the state £ k . The data points correspond to sev- 
eral measurements binned in 1 ms intervals and error bars 
representing ± the standard deviation. The blue/red/green 
curves correspond to measurements at different offset fields 
with quadratic Zeeman periods of ir/q « 23/11.7/5.7 ms re- 
spectively. The lines correspond to the analytic solution ||1J, 
with parameters q calculated from the magnetic offset field 
and k obtained from fitting to the first 10 ms of data. Ad- 
ditional fit parameters account for detection errors (y-offset) 
and evolution during time of flight (x-offset). 




calculated |k| 

FIG. 4: Fit values of the |£o| 2 oscillation amplitude parame- 
ter k of Eq. Mat versus the theoretical expectation for mea- 
surements with 25 ... 30 % condensate fraction (squares) and 
with 75 ... 80 % condensate fraction (triangles). The lines cor- 
respond to linear fits excluding the highest k data points. 



tal results to our theoretical prediction, we plot the fit 
value (fitting the analytic solution to the first 10 ms of 
data) versus the calculated value (using the scatte ring 
lengths a F =2 = 100.4(1) o and a F=0 = 101.8(2) a Hf 
and the Thomas- Fermi approximation to obtain (n)) in 
Fig. 01 The fitted values indicate a lower effective in- 
teraction parameter of 0.2 to 0.4 times the calculated 
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value and surprisingly show a clear temperature depen- 
dence. Possible explanations include the breakdown of 
the Thomas-Fermi approximation in the tightly confined 
vertical direction, trap anharmonicities or small scale do- 
main formation and require further investigations. 

In conclusion we have analyzed the interplay of exter- 
nal magnetic fields and interatomic interactions on the 
system evolution, giving analytic expressions in compar- 
ison to experimental data. The experimental data at 
finite temperature shows various dephasing/decoherence 
mechanisms, which we were able to separate using the 



theoretical model and adapted experimental sequences. 
These investigations pave the way to use spinor conden- 
sates for general studies of dephasing and decoherence 
phenomena in multi-component systems. 
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